Double strand breaks in DNA resulting from double-electron-emission events 
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A mechanism of double strand breaking (DSB) in DNA due to the action of two electrons is con- 
sidered. These are the electrons produced in the vicinity of DNA molecules due to ionization of water 
molecules with a consecutive emission of two electrons, making such a mechanism possible. This 
effect qualitatively solves a puzzle of large yields of DSBs following irradiation of DNA molecules. 
The transport of secondary electrons, including the additional electrons, is studied in relation to the 
assessment of radiation damage due to incident ions. This work is a stage in the inclusion of Auger 
mechanism and like effects into the multiscale approach to ion-beam cancer therapy. 

PACS numbers: 87.53.-j, 61.80.-x, 34.50.Gb, 36.40.Cg 



The analysis and assessment of radiation damage are 
important in relation to a wide range of applications from 
cancer therapy to radiation protection of humans and 
electronics. Studies of biodamage due to irradiation with 
ion beams, centered on the analysis of pathways of DNA 
damage, are discussed in a large number of papers ob- 
served in reviews 0, Q ; the multiscale approach to the 
physics of ion-beam cancer therapy [3|, |J] is also among 
these studies. A double strand break (DSB) in DNA is 
the most pernicious kind of damage that happens as a 
result of radiation impact. It is defined as the simultane- 
ous breaking of the two strands of a DNA molecule within 
the distance of 10 base pairs along the helix, which corre- 
sponds to a single convolution of a DNA molecule. This 
type of lesion is emphasized because of the difficulties 
of its repair and thus close connection to cell lethality. 
What processes lead to this lesion? 

If the primary ionizing projectiles are ions, a substan- 
tial fraction of damage is done by secondary electrons, 
formed in the process of ionization of the medium. A 
number of different pathways of damage due to these 
electrons were considered in Refs. [H-[(|- The direct mea- 
surements of DNA damage due to incident electrons were 
presented in, e.g., Ref. [7[, which brought to light the pos- 
sibility that low energy electrons were important agents 
of DNA damage. Since then, the mechanism of a SSB due 
to the action of a single electron, related to the formation 
of a transient negative ion (TNI) as a part of the pro- 
cess of dissociative electron attachment (DEA) has been 
widely discussed in the literature, e.g., Refs. @ with 
emphasis on low (under ionization threshold) energy of 
the incident electrons. 

Unexpectedly high yields of DSBs compared with SSBs 
in Ref. 01 lead to a hypothesis that DSBs can be caused 
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by a single electron. This logic has been widely accepted 
and used by the present authors to calculate the yields of 
DSBs due to ions @. However, the mechanism of DSBs 
due to low-energy electrons is still quantitatively unclear 
despite qualitative arguments, suggesting that the breaks 
in the second strand are due to the action of debris gen- 
erated by the first SSB ■ In Ref. Q , the production of 
DSBs by two separate electrons was also considered, but 
that analysis was then shelved, since the number density 
of secondary electrons due to primary ionization with an 
ion was not nearly enough for this effect to be consid- 
erable in comparison with DSBs due to single electron 
action. 

In this work, we argue that additional electrons emit- 
ted in the vicinity of a DNA molecule as a result of the 
Auger-like mechanism augment the above effect and thus 
constitute a mechanism for a DSB. We explore the prob- 
ability of two electrons, produced in the vicinity of a 
DNA molecule, be incident on a single convolution of 
this molecule. These additional electrons emerge as a re- 
sult of double ionization events such as the Auger effect in 
single molecules, e.g., Refs. for water, or to the ef- 

fect of intermolecular coulombic decay (ICD), studied in 
Refs. HUH 

for water clusters. The actual mechanism 
of double ionization is not important for our analysis. 
We will use a generic term "double-ionization-events" to 
describe all relevant events leading to a production of 
additional electrons, and we will refer to the additional 
electron as an Auger electron, even if it originates from 
ICD. 

1. Auger electrons are produced as a result of non- 
radiative relaxation of holes produced by primary and 
secondary ionization. They may emerge consequent to 
the ionization of water or other molecules or clusters of 
the medium. If a secondary electron ionizes a molecule 
by kicking out one of its inner-shell electrons, a hole is 
formed on this molecule. If this hole then relaxes via a 
non-radiative channel (an Auger electron is emitted from 
the same molecule or an adjacent molecule in a cluster) 
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FIG. 1: The scheme of ionization with the Auger mechanism 
applied to ICD; e2 and the Auger electron emerge from dif- 
ferent molecules. The geometry of the problem; ionization 
happens at the origin; the cylinder represents a DNA convo- 
lution. 



the total number of electrons, emerging from a sub-nm lo- 
cality, is equal to three, the ionizing electron, the released 
electron, and the Auger electron, as shown in Fig. [TJ If a 
double-ionization event occurs on a DNA molecule, the 
probability of damage, such as a DSB, comprises the se- 
quence of these processes with a corresponding dynam- 
ics of the DNA molecule and possible further interaction 
with these three electrons. If a water molecule (cluster) 
of the medium located in the vicinity of a DNA molecule 
is a host of a double-ionization event, then these three 
electrons can independently diffuse and stumble onto a 
DNA convolution and produce two SSBs, possibly yield- 
ing a DSB. Let us consider the latter scenario in more 
detail. 

It is interesting to notice that the energies of each of 
the three electrons, emerging from a water molecule or 
cluster, are likely to be below 15 eV. An average ionizing 
electron has the energy around 45-50 eV, corresponding 
to the average energy of electrons produced by an io n 1191 . 
After this ionization event, it loses about 30 eV [la]- 1 
Thus, the ionizing electron is left with less than 15 eV 
while the other two will have even smaller energies [l6| . 
These low-energy electrons are likely to be engaged into 
the DEA channel leading to SSBs. How does this scenario 
make a difference in the DSBs production? 

The number density of secondary electrons, produced 
on the ion's path, reduces rather steeply with the in- 
creasing distance from the path, so that the probability 
of two electrons incident on a single convolution of a DNA 



1 In order to kick out an electron from la\ MO of H2O molecule, 
the secondary electron has to have the energy above 540 eV. Such 
electrons ere very rare. Ionization of other state cannot produce 
Auger electron because of the lack of energy. 



molecule, located a few nm from the path is very small. 
However, when ionization of a water molecule (cluster) 
producing an Auger electron occurs at a distance from 
the path, three electrons emerging from the same spot 
substantially boost the local number density of electrons 
and, hence, the probability of a nearby DNA convolu- 
tion to be hit with two electrons. A similar process may 
take place if the primary projectile is a photon. In that 
case, two (instead of three) electrons are produced in a 
locality and are capable of producing a DSB in a DNA 
molecule. If incident photons ionize the lai state, then a 
cascade of Auger electrons may follow, rapidly increasing 
the number density of electrons. 

In order to give a quantitative example of the effect 
caused by the Auger mechanism, let us consider an event 
of ionization with Auger emission caused by a secondary 
electron, produced by an incident carbon ion. The goal 
of this example is to calculate the transport of three elec- 
trons emerging from the ionization locality, which we will 
treat as a point, to a nearby DNA convolution. 

Let an event of double ionization of a water molecule 
(cluster) happen at the origin. We consider a three- 
dimensional random walk of electrons from this point. 
We represent a single DNA convolution with a cylinder of 
radius a = 1.15 nm and a length of 3.4 nm. For simplicity 
and definiteness, we situate this cylinder symmetrically 
with respect to the y-axis, with the cylinder's axis par- 
allel to the x-axis and lying in the xy plane, as shown in 
Fig.ffl The distance between the origin and the cylinder's 
axis, p, exceeds a so that the Auger electron is emitted 
outside the convolution. 

If one electron is emitted from the origin at t = 0, 
according to Ref. [2(J, its rate, dpi/dt, of passing through 
the patch, cL4, located at a distance r from the origin, is 
given by the expression 



d ^= d A.Dn T dP ^ r) 



dt 



Or 



(1) 



where D = vl/6 is the diffusion coefficient, v is the speed 
of the electron, n r is a unit vector in the radial direction, 
and 



P(t,r) 



p_f\ xp (_iri) 



(2) 



is the probability density to observe a randomly walking 
electron at a time t and a distance r from the origin. 
Eq. ([l} should be integrated over both the time and dA, 
in order to calculate the probability for the electron to 
encounter the cylinder. 

The time dependence in Eqs. (fl] [2]) can be translated 
into the dependence on number of steps, k, using vt = kl, 
where I is the elastic mean free path of electrons in the 
medium. Then, we rewrite Eq. (fl}, substituting @, and 
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FIG. 2: The probability for two electrons to pass through a 
single convolution of DNA (solid line) compared to that of 
one electron (dashed line). 



switching from variable t to k as 



dpi (r, k) = dkdA ■ n, 



r 

2k 



2nkf 



3/2 



exp 



3r 2 
2kf 



(3) 



and integrate it over k (from the minimal number of 
steps necessary to reach the surface of the cylinder to 
infinity) and dA over the surface of the cylinder. As 
was noticed above, these electrons are very near or be- 
low the ionization threshold and their interactions with 
water molecules will be, by and large, elastic. Therefore, 
we integrate Eq. [3] without including attenuation due to 
inelastic collisions and using the value for the mean free 
path, 1 = 0.15 nm dJ. 

The results of this integration of are presented in Fig. [5] 
The probability decreases with increasing r. This proba- 
bility is proportional to that of producing a strand break 
or a different type of lesion on the length of a given DNA 
convolution. The corresponding coefficient is the proba- 
bility of inducing an SSB in a single electron action, T, 
remains largely unknown and its calculations and mea- 
surements remain a task for atomic physicists and quan- 
tum chemists. 

The second and third electrons from the three, pro- 
duced in the double ionization event, are independent 
from the first one and, hence, the probability of two elec- 
trons to impact the cylinder, p 2 , is given by 



p 2 (r) = Zp\(r) 



(4) 



The dependence of this probability on the distance of 
the cylinder axis from the origin is shown in Fig. [2] It 
is remarkable, that for r < 2nm the probability of the 
impact of two electrons is comparable to that of one. 
Of course, Fig. [5] includes neither the value of T, nor 
the probability of an ionization event with the following 



Auger emission, but still, a substantial quantity of this 
fluence makes the Auger-mechanism influence on DSB 
yield viable. 

2. Auger electrons emitted from the ion's path as a result 
of primary ionization propagate similarly to secondary 
electrons, so the correction to the number of agents in- 
teracting with a DNA molecule situated at a distance 
from the path can be made after the calculation of the 
probability of Auger electron emission during the pri- 
mary ionization. Our goal in this section is to consider 
the Auger electrons emitted as a result of the ionization 
of water molecule clusters by secondary electrons. 

The scenario is as follows: Secondary electrons orig- 
inate on the ion's path and diffuse (through a random 
walk) away from the path. Their interactions with water 
molecules are primarily elastic until they ionize a wa- 
ter molecule (cluster). This happens about once per ev- 
ery 30 elastic collisions at relevant energies of secondary 
electrons. A typical secondary electron ionizes one wa- 
ter molecule before it becomes thermalized and finally 
bound. More rare energetic 5-electrons are not consid- 
ered in this discussion. For simplicity, let us assume that 
ionization happens on a certain step in the random walk 
of a 45-eV electron. At that time, this electron is most 
likely situated at a distance from the path, given by the 
expression, 



pP(k,r)d 3 r 



(5) 



where k — 30 and p is the radial distance from the 
path. The calculations, with the above parameters for 
secondary electron propagation in water, give p approx- 
imately equal to 1 nm. This means that the locus of 
first ionization events is a cylinder with a radius of 1 nm 
on the axis of the ion path. Then, we can calculate the 
fluence at a distance p from the path due to Auger elec- 
trons emitted from the surface of this cylinder. It is given 
by the integration of Eq. ^ normalized by dA over the 
surface of the cylinder: 



Pa(p) 





r 


I 


~ ■ n P 




r 
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where r 2 = p 2 



+ ( 2 + p 2 — 2pp cos (j>, C is the coordinate 
along the path, cj) is the azimuthal angle, and n p is a unit 
vector in the radial direction. The rest of the geometry 
is shown in Fig. [3] along with the results for pa(p)- The 
fluence per one electron from the cylinder is compared to 
that of electron coming from the path. Starting from a 
small distance from the cylinder, the fluences are equal, 
which means that outside some domain, the total fluence 
is equal to the fluence of secondary electrons emitted on 
the path multiplied by (1+2?/;), where ip is the probability 
of the double-ionization event on impact and the factor 
of two is due to an extra electron emitted in this process. 
3. The quintessence of this paper is that due to the Auger 
mechanism, events in which two or more electrons inter- 
act with a single DNA convolution are not rare. Double 
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FIG. 3: The dependence of the fluence of an Auger electron 
from ionization by a secondary electron on the distance from 
the path (solid line) compared to that of a single electron 
emitted from the path (dashed line). 

strand breaks or other types of complex damage may 
result from these interactions. Thus, this effect quali- 
tatively provides a solution to the puzzle of high yields 
of DSBs by finding sources of high local electron den- 
sity near sites of damage. Since the Auger emission may 
result from the ionization of a water molecule by an inci- 
dent photon, the above two-electron mechanism can also 



contribute to the DSB yield consequent to an irradiation 
of tissue with photons. 

We considered the major concepts related to the trans- 
port of Auger electrons and made first steps to includ- 
ing the Auger electrons in the multiscale approach to 
the calculation of radiation damage by ions. Besides 
the consideration of a two-electron mechanism for DSBs, 
our findings also correct the fluence of secondary elec- 
trons through the DNA convolution due to the Auger 
mechanism. The inclusion of this effect into the mul- 
tiscale approach is an important step, since, in general, 
the Auger mechanism plays a significant role in radiation 
damage [22| . 

Future research in this route will be based on the 
calculations and measurements of the cross sections of 
ionization of water molecules and clusters with a conse- 
quent emission of Auger electrons and the cross sections 
of DSBs due to two incident electrons. 

The calculation of the fluence of Auger electrons pre- 
sented in the previous section gives a framework for the 
calculation of the transport of free radicals formed due 
to an ion's traverse through a medium. This calculation, 
without the use of Monte Carlo simulations has long been 
desired. Finally, the analysis presented in this paper can 
be extended to the calculation of complex damage as in 
Ref. Q with the inclusion of the Auger mechanism. 
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